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1 Introduction 

The standard model (SM) of particle physics has been successful in explaining a wide variety of 
data. In spite of this, the SM is incomplete. For example, it possesses a divergence in the Higgs 
sector [1 J and has no cold dark matter (DM) candidate [2j. Many models of physics beyond the 
SM (BSM) have been proposed in order to address these problems. 

DM particles, if produced in proton-proton collisions at the CERN Large Hadron Collider 
(LHC), would escape detection and result in a significant transverse momentum (pj) imbal- 
ance in the detector. Additionally, cascade decays of heavy colored particles to final states with 
a high multiplicity of energetic jets and T leptons appear very naturally in many BSM physics 
scenarios. Hence, events with multiple T lepton candidates, large jet multiplicity, and signif- 
icant transverse momentum imbalance, represent a distinct signature of new physics. In this 
paper, focus is placed on final states with hadronically decaying t leptons. In what follows, the 
visible part of a hadronically decaying t lepton will be referred to as T^ 

In certain models of supersymmetry (SUSY), the lightest supersymmetric particle (LSP) is a 
candidate for DM. It has been appreciated for some time that the DM relic density may be 
sensitive to coannihilation processes involving the LSP and the next-to-lightest supersymmetric 
particle (NLSP). Coannihilation is characterized by a mass difference (AM) between the NLSP 
and the LSP of approximately 5-15 GeV [3-6] . This small mass difference would be necessary 
to allow the NLSP to coannihilate with the LSP in the early universe, leading to the dark matter 
abundance that is currently observed [7]. If the supersymmetric partner of the T lepton, the 
stau (f ), is the NLSP, and if the f decays primarily to a x lepton and the LSP, small values of 
AM would lead to final states with low-energy t leptons (px ~ AM) (H. Decays of colored 
SUSY particles can produce the f via chargino (^ ± ) or neutralino (x°) intermediate states (e.g., 
^2 — > tt — > tt^j 1 ), resulting in final states with at least one T^. 

We present a search for BSM particles in events with exactly one Th lepton and jets (single-Th 
final state), and in events with jets and two or more leptons (multiple-Th final state). These 
two topologies provide complementary sensitivity to models with a wide range of AM values. 
For example, in the case of very small values of AM (~5GeV), the low-energy Th cannot be 
effectively detected and the search for new physics in the single-Th final state has better sen- 
sitivity. The analysis is performed using proton-proton collision data at i/s = 7 TeV collected 
with the Compact Muon Solenoid (CMS) detector [9] at the LHC in 2011. The data sample 
corresponds to an integrated luminosity of 4.98 ± 0.11 fb _1 . The search is characterized by 
methods that determine the backgrounds directly from data, to reduce the reliance on simu- 
lation. To illustrate the sensitivity of this search to BSM processes, the constrained minimal 
supersymmetric extension of the standard model, or minimal supergravity is chosen as the 
benchmark (3][l0l[TTl; we denote this benchmark as "CMSSM". An interpretation of the results 
in the context of simplified model spectra (SMS) |l2[[l3l is also presented. 

2 The CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid, of 6 m inner diameter, 
providing a magnetic field of 3.8 T. Within the field volume are a silicon pixel and strip tracker, 
a crystal electromagnetic calorimeter (ECAL), which includes a silicon sensor preshower de- 
tector in front of the ECAL endcaps, and a brass-scintillator hadron calorimeter. Muons are 
measured in gas-ionization detectors embedded in the steel return yoke. In addition to the 
barrel and endcap detectors, CMS has extensive forward calorimetry. 
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3 Object reconstruction and identification 



The inner tracker measures charged particles within \rj\ < 2.5 and provides an impact parame- 
ter resolution of about 15 /mi and a pj resolution of about 1.5% for 100 GeV particles. Collision 
events are selected with a first-level trigger based on fast electronics, and a higher-level trigger 
that runs a version of the offline reconstruction program optimized for speed. 

The CMS experiment uses a right-handed coordinate system, with the origin at the nominal 
interaction point, the X axis pointing to the center of the LHC ring, the y axis pointing up 
(perpendicular to the plane of the LHC ring), and the z axis along the counterclockwise beam 
direction. The polar angle 6 is measured from the positive z axis and the azimuthal angle in the 
x-y plane. The pseudorapidity is given by rj = — ln[tan(0/2)]. 

3 Object reconstruction and identification 

Jets in the detector are reconstructed using particle-flow (PF) objects [14J. In the PF approach, 
information from all subdetectors is combined to reconstruct and identify final-state particles 
(muons, electrons, photons, and charged and neutral hadrons) produced in the collision. The 
anti-fcx clustering algorithm [15] with a distance parameter R = 0.5 is used for jet clustering. 
Jets are required to satisfy criteria designed to identify anomalous behavior in the calorimeters, 
and to be well separated from any identified x lepton. 

Validation and efficiency studies are performed utilizing events with a Th lepton and a light- 
lepton t, with I representing an electron (e) or muon Muons are reconstructed using the 
tracker and muon chambers. Selection requirements based on the minimum number of hits 
in the silicon tracker, pixel detector, and muon chambers are applied to suppress muon back- 
grounds from decays-in-flight of pions or kaons |[T6ll . Electrons are reconstructed by combining 
tracks with ECAL clusters. Requirements are imposed to distinguish between prompt and non- 
prompt electrons, where the latter can arise from charged pion decay or photon conversion [fl~7|. 
The light-lepton candidates are required to satisfy both track and ECAL isolation requirements. 
The track isolation variable is defined as the sum of the pi of the tracks, as measured by the 
tracking system, within an isolation cone of radius AR = yj (A//) 2 + (Acp) 2 = 0.4 centered on 
the light-lepton track. The ECAL isolation variable is based on the amount of energy deposited 
in the ECAL within the same isolation cone. In both cases the contribution from the light-lepton 
candidate is removed from the sum. 

Reconstruction of hadronically decaying x leptons is performed using the hadron-plus-strips 
(HPS) algorithm [18|, designed to optimize the performance of Th reconstruction by consider- 
ing specific Th decay modes. To suppress backgrounds in which light-quark or gluon jets mimic 
hadronic T decays, a Th candidate is required to be spatially isolated from other energy deposits 
in the calorimeter. Charged hadrons and photons not considered in the reconstruction of the Th 
decay mode are used to calculate the isolation. Additionally, candidates are required to be 
distinguished from electrons and muons in the event. In this analysis, two HPS isolation defi- 
nitions are used. The Th isolation definition used for single-Th final states rejects a Th candidate 
if one or more charged hadrons with pj > 1.0 GeV or one or more photons with transverse 
energy Ej > 1.5 GeV is found within an isolation cone of radius AR = 0.5. The Th isolation def- 
inition used for multiple-Th final states rejects a Th candidate if one or more charged hadrons 
with pi > 1.5 GeV or one or more photons with transverse energy Ej > 2.0 GeV is found within 
an isolation cone of radius AR = 0.3. The isolation criteria used for the multiple-Th final state 
increases the signal-to-background ratio while reducing the rate of Th misidentification. This 
affects the yield of events with light-quark or gluon jets that are misidentified as Th leptons, 
which depends on the square of the misidentification rate. 



3 



The missing transverse momentum #fr is defined as: 



(i) 



where the sum runs over all the jets in the event with pj > 30 GeV. The Iflj vector is the 

negative of the vector sum in Eq. The observable Hj = YIPj is used to estimate the 
overall energy scale of the event. For the single-Th final state, Hj is calculated using jets with 
pj > 50 GeV and will be referred to as H|°. For the multiple-^ final state, Hj is calculated 
using jets with pj > 30 GeV and will be referred to as Hj°. The sum #fr + HS will be referred 
to as M e ff . The use of a lower pj threshold for the jets in the multiple-Th final state increases the 
efficiency of signal events without significantly increasing the background. 



4 Signal and background samples 

The major sources of SM background are top-quark pair (tt) events and events with a W or 
Z boson accompanied by jets. Both tt and W + jets events can have genuine Th leptons; large 
genuine Iflj from W boson decays; and jets that can be misidentified as a Th. Similarly, Z + jets 
events with Z(— > vv) and with one or more jets misidentified as a Th lepton provide a source 
of background. Z + jets events with Z(— > vv) present a background because of the genuine 
Th leptons and the genuine Jflj from the neutrinos in the Th decay. QCD multijet events can 
become a background when a mismeasured jet gives rise to large and jets are misidentified 
as a Th lepton. 

Data are compared with predictions obtained from samples of Monte Carlo (MC) simulated 
events. Signal and background MC samples are produced with the pythia 6.4.22 [ 19 1 and 
MadGraph [20 1 generators using the Z2 tune EU and the NLO CTEQ6L1 parton distribution 
function (PDF) set |22"|. The t lepton decays are simulated with the TAUOLA |23| program. 
The generated events are processed with a detailed simulation of the CMS apparatus using the 
Geant4 package [24]. The MC yields are normalized to the integrated luminosity of the data 
using next-to-leading order (NLO) cross sections ff25l430| . For the 2011 LHC running condi- 
tions, the mean number of interactions in a single beam crossing is ~10. The effect of multiple 
interactions per bunch crossing (pileup) is taken into account by superimposing MC minimum- 
bias events so that the probability distribution for overlapping pp collisions in the simulation 
matches the measured distribution. 



5 Event selection 

Events for both the single- and multiple-Th final states are selected using a trigger that requires 
> 150 GeV. This trigger allows us to maintain sensitivity in regions where the pj value 
of the Th is small (pj ~ 15 GeV). This trigger efficiency, for an offline selection requirement of 
Iflj > 250 GeV, is 98.9%. For the Th efficiency and validation studies, samples are chosen using 
triggers that require the presence of both a Th candidate and a muon. 

The Th candidates must satisfy pj > 15 GeV and \rj\ < 2.1. For the single-Th final state we 
require that no additional light leptons be present in the event. This requirement suppresses 
background from tt, W + jets, and Z + jets events. For the multiple-Th final state there is no 
requirement placed on the number of light leptons. 

For the single-Th final state, we define a baseline event selection H|° > 350 GeV and > 
250 GeV. The sample obtained with the baseline selection is used to validate the background 
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6 Background estimate 



predictions. The signal region (SR) for the single-ih final state is defined by H|° > 600 GeV and 
# T >400GeV. 

For the multiple-Th final state, the SR is defined by #fr > 250 GeV and by the requirement that 
there be at least two jets with pj > 100 GeV and \rj \ < 3.0. QCD multijet events are rejected by 
requiring the azimuthal difference A<^(j2, Vlj) between the second leading jet in pj and the Ifij 
vector to satisfy | ^(p{ji, Vlj) | > 0.5. Finally events are required to contain at least one t^t^ pair 
separated by AR(T h ;/ r h y) > 0.3. 



6 Background estimate 

The background contributions are categorized differently for the single- and multiple-Th final 
states. For the single-Th final state, the background contributions are divided into events con- 
taining a genuine and events where a jet is misidentified as a t^. For the multiple- final 
state, the main background contribution arises from misidentified Th leptons. We identify the 
different sources of background individually using dedicated data control regions (CR). 

6.1 Estimate of backgrounds in the single-r h final state 

In the single-Th final state, the largest background contribution comes from W + jets events that 
contain a genuine Th lepton. The other significant contribution arises from QCD multijet events 
in which a jet is misidentified as a Th- The W + jets background contribution is estimated using 
a sample of W + jets events with W — > }iv. The QCD multijet background is determined by 
selecting a QCD-dominated CR and evaluating the Th misidentification rate. 

6.1 .1 Estimate of the W + jets background in the single-r h final state 

To evaluate the W + jets background, we exploit the similarity between W decays to a muon 
and to a tau lepton and select a sample of W + jets events with W(— > ]iv). This sample will be 
referred to as the muon control sample. To select the muon control sample, events are required 
to contain exactly one muon and no reconstructed Th or electron. To emulate the Th acceptance, 
the muon is required to satisfy \t]\ < 2.1. The yields in the muon control sample are corrected 
for muon reconstruction (e^ eco ) and isolation efficiency (e^ ). The muon reconstruction effi- 
ciency is derived from data using a sample of Z + jets events and parameterized as a function 
of pj and rj. The muon isolation criteria help to distinguish between muons from the decay of 
the W boson and muons from semileptonic decays of c and b quarks. The isolation efficiency 
is parameterized as a function of the separation from the nearest jet and the momentum of the 
jet. A correction factor (P™) is applied to the muons in the muon control sample to account for 
muons that do not come from a T-lepton decay. This correction factor depends on the pj of the 
muon and the tflj value in the event and is derived from a simulated sample of W + jets events. 

As the muons in the muon control sample are selected to mimic a Th, a correction is applied to 
emulate the probability to reconstruct and identify a Th lepton. The reconstruction and iden- 
tification efficiency e^ eco is parameterized as a function of the pj of the Th candidate and as a 
function of the total number N of charged particles and photons in the isolation cone [Fig. [l^] . 
Corrections are also applied to account for the hadronic branching fraction (y£ f ( hadr )) Q f a t 
lepton. Except for the yk f ( hadr ) the values of the correction factors differ in each event. The 
corrections are combined to define an overall event weight, defined as: 

PW x e T x f* f(hadr) 

rcorr _ f* ■> iy>. 

J event — £reco x £ i so • V-) 



6.1 Estimate of backgrounds in the single-r h final state 
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Figure 1: (a) Dependence of the Th reconstruction efficiency e r ^ co on the number of additional 
particles N in the isolation cone in bins of Th lepton pj for the single-Th final state, where N is 
the total number of the photons and charged hadrons in the isolation cone, and (b) dependence 
of Pj gen on the Th response. Both distributions are derived from a simulated sample of W — > 
rv) + jets events. 

A Th response template is derived from simulated events. The response template is given by the 
ratio of the reconstructed energy of the Th to the true generator-level energy. The Th response 
depends on the transverse momentum of the generated t lepton [Fig. [lj?] and on the number 
of reconstructed primary vertices in the event. The muon pj spectrum is smeared as a function 
of pj and the number of primary vertices to mimic the pj distribution of the Tjj. 

Fully simulated W + jets events are used to verify the procedure. Figure [2] shows the H|° 
and Vlj distributions from simulated W + jets events for the single- Th final state. These events 
satisfy the baseline selection described in Section [5] The reconstructed Th is required to match 
a hadronically decaying generated tau lepton, to ensure that only the genuine tau background 
is addressed in this check. The event yield and distributions are compared with the prediction 
from the muon control sample and agree within statistical uncertainties, thus verifying the 
closure of the method in MC simulation. Hence, the predicted H T ° and distributions from 
the muon control sample can be taken to describe a Th sample within statistical uncertainties. 

Table 1: The selected and predicted background contributions for simulated events with a gen- 
uine Th passing the baseline and signal selection in the single- Th final state. The reconstructed 
Th is required to match the visible part of the generated, hadronically decaying T-lepton. The 
predictions are derived from the muon control sample. 



L = 4.98fb~ 



Baseline Selection Signal Selection 

Selected Predicted Selected Predicted 

W(-> lv) +jets 452 ± 30 441 ± 21 28.9 ± 7.5 34.9 ± 5.9 

tt 60.6 ±3.7 63.2 ±2.1 1.6 ±0.6 2.9 ±0.4 

Z(-»#) + jets 10.9 ±2.1 8.4 ±1.3 0.8 ±0.6 0.4 ±0.3 

W+W~ 15.1 ±1.6 14.4 ±1.1 0.5 ±0.3 1.3 ±0.3 

Sum 539 ±30 527 ±21 31.8 ±7.5 39.5 ±5.9 



The muon control sample consists primarily of W + jets events, but also contains tt events in 
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6 Background estimate 




Figure 2: Distributions of (a) Hj and (b) I/It for the genuine Th estimate in simulated W + jets 
events for the single-Th final state. The black triangles show the results for events that satisfy the 
baseline selection and that contain a reconstructed Th matched to the visible part of a generated, 
hadronically decaying t lepton. The filled green areas show the prediction obtained from the 
muon control sample. The hatched areas are the total uncertainty on the prediction. 

which one W boson boson decays into a muon while the other W boson decays either into 
an unidentified T lepton or into a light lepton that is not reconstructed. Any isolated muons 
produced through the decay of b or c quarks can also contribute to the muon control sample. 
SM processes containing a Z boson or two W bosons can also contribute to the muon control 
sample if one of the two decay muons is not reconstructed. 

The true event yields of each process as determined from simulation are summarized in Table [T] 
for the baseline and SR selections. For both selections the number of predicted events with a 
genuine Th lepton is seen to agree with the true number of events. The value of e^ eco that is 
used to calculate the predicted rate is measured in a sample of W + jets and is different from 
the value that would be measured in a sample of tt events. This leads to an overestimation of 
the tt contribution. A systematic uncertainty is assigned to account for this overestimation. 

6.1 .2 Estimate of the QCD multiet background in the single-r h final state 

To estimate the background where a jet is misidentified as a Th lepton, a QCD-dominated con- 
trol sample is obtained by selecting events with Hj° > 350 GeV and 40 < Wt < 60 GeV. The 
control sample is selected using a prescaled Hj trigger with criteria that lead to a sample where 
about 99% of the events arise from QCD multijet production. The probability for a jet to be 
misidentified as a lepton is measured by determining the fraction of jets from the single-Th 
control sample that pass the Th identification criteria. Jets considered in the calculation of the 
misidentification rate satisfy the requirements pj > 5 GeV and \tj\ < 2.5. The misidentification 
rates /; for each jet i depend on t] and px arid are used to determine an overall weight, which is 
applied to each event. The event weights are defined as: 

<ent = l-fT(l ( 3 ) 
i 

where n is the number of jets. The measured misidentification rates shown in Fig. [3](a) are 
applied to data events in the region with H T ° > 350 GeV and with H T ° > 600 for two regions 



6.2 Estimate of backgrounds in the multiple-T h final state 
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Figure 3: (a) The rate of jet misidentification as a lepton in simulation (triangular symbols) 
and data (circular symbols) as a function of p J j* for events with H5° > 350 GeV and 40 < yij < 
60 GeV; (b) The Vlj distribution estimated in simulated events with H|° > 350 GeV, where the 
triangular symbols represent events that pass the baseline selection, the filled blue area shows 
the predicted events, and the hatched area shows the total uncertainty on the prediction. These 
distributions correspond to the single-Th final state. 

of fir- 60 < Vl T < 80 GeV and 80 < # T < 100 GeV. These four regions are dominated by 
QCD multijet events. The results for data and simulation, as well as the predicted fraction 
of QCD multijet events, are shown in Table [2] The ratio of selected events over predicted 
events is statistically compatible with one and stable over the range of Iflj. Figure |3jb) shows 
the Iflj distributions of predicted and selected events for simulated QCD multijet events with 
Hj° > 350 GeV. The two distributions agree over the whole range of Wt ■ 

Table 2: The percentage of QCD multijet events in the Vlj binned samples for different QCD 
multijet dominated regions in the single-Th final state. 







Vlj [GeV] 




Hf > 350 GeV 


60-80 


80-100 


> 250 


QCD fraction 


97% 


93% 


6% 


selected /predicted (sim) 


0.98 ± 0.06 


0.96 ± 0.07 


1.24 ±0.28 


selected/ predicted (data) 


1.01 ±0.08 


0.88 ±0.13 








Wt [GeV] 




H|° > 600 GeV 






>400 


QCD fraction 


96% 


93% 


17% 


selected /predicted (sim) 


0.94 ± 0.09 


0.85 ± 0.09 


2.43 ± 1.45 


selected/ predicted (data) 


1.14 ±0.26 


0.97 ±0.37 





6.2 Estimate of backgrounds in the multiple-Th final state 

The estimate of the SM background contributions to the SR sample for multiple-Th events is 
based on the number of observed events in CRs. The events in each CR are selected with 
similar selection requirements to those used in the SR, but are enriched with events from the 
background process in question. Correction factors and selection efficiencies are measured in 
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those CRs and used to extrapolate to the SR. We use the observed jet multiplicity in each CR 
along with the measured rate at which a jet is misidentified as a Th to calculate the yield in the 
SR. The following equation is used to estimate each background contribution B: 

Nf = N^[k tt V(0) + cc rj V(l) + ec jS P(2)] f (4) 

where N| R is the predicted rate in the SR, Ng R is the observed number of events in the CR, 
and d X y is the correction factor for acceptance and efficiency for events in the CR with true 
physics objects "x" and "y" . Here the physics object can be a Th or a quark or gluon jet. Since 
the dominant SM backgrounds contribute to the SR when jets are misidentified as Th lepton, 
the background estimation strategy outlined in Eq. relies on the determination of the event 
probability Vim) for at least "m" jets to be misidentified as a Th, where Vim) is the product of 
three factors: (i) the probability P(N) for an event to contain N jets, (ii) the number of possible 
ways for exactly n jets to pass the Th identification criteria given N possible jets C(N,n) = 
N\/n\(N — n)!, and (iii) the probability / for a single jet to be misidentified as a T^. The Vim) 
terms are given by: 

co N 

Vim) = £ P(N) E C(N,ft)r(l-/) N -". (5) 

N=m n=m 

Equation (|5j would be identical to Eq. (|3j if used in the case of the single-Th final state. Equa- 
tion (jij) is used to estimate the tt, W + jets, and Z + jets background contributions to the SR. The 
P(N) terms are determined from data using the jet multiplicity distribution in each CR, while 
the / terms are measured for each background process by determining the fraction of jets in 
each CR that pass the Th identification criteria. Since the QCD multijet contribution to the SR 
for the multiple-Th final state is negligible according to simulation, a data-to-MC scale factor 
is used to correct the QCD multijet prediction from simulation. In the sections that follow, the 
selections used to define high purity CRs are outlined and the correction factors a xy used in 
Eq. Q are defined. The fraction of events with two Th leptons is denoted A TT , the fraction with 
one Th lepton and one jet misidentified as a Th lepton is denoted A T j, and the fraction with two 
jets misidentified as Th leptons is denoted Au. 

6.2.1 Estimate of the tt event background to the multiple-r h final state 

To estimate the contribution of tt events to the multiple-Th SR, a CR is selected by removing 
the Th isolation requirement and by requiring the presence of at least two b-quark jets (b jets), 
identified using the track-counting-high-efficiency (TCHE) algorithm at the medium working 
point II3T1 . Because QCD multijet, W + jets, Z(— >■ tt) + jets and Z(— > vv) + jets events are 
unlikely to contain two b jets, this requirement provides a sample in which about 99% of the 
events are tt events, according to simulation. Figure [IJa) shows the px distribution of Th leptons 
in the tt CR for data and simulation. 

According to simulation, the fraction of events in the tt control sample that contains one gen- 
uine Th is A T j = 0.166 ± 0.011, while the fraction without a genuine Th is Au = 0.834 ± 0.025. 
The genuine ThTh contribution is negligible (A TT ~ 0) according to simulation. Incomplete 
knowledge of the genuine ThTh contribution is included as a source of systematic uncertainty 
in the tt background prediction. Therefore, cc X j in Eq. (jijl is given by A T ye^°/ P(2 b jets), where 
e!^° is the probability for a Th lepton to pass the isolation requirement, while ocu is given by 
Ayy/P(2bjets). The probability P(2 b jets) to identify two or more b jets is determined by the b 



6.2 Estimate of backgrounds in the multiple-T h final state 
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jet identification efficiency factor pl|. The number of tt events in the SR is calculated as: 

The probability for a jet in a tt event to be misidentified as a Th lepton has an average measured 
value of / = 0.022 ± 0.004. Cross checks are made to validate the use of the b-jet identification 
efficiency as measured in Ref. f3lf for this analysis. The estimated tf contribution in the SR is 
determined to be = 2.03 ± 0.36. 




200 400 600 800 1000 400 600 800 1000 



Figure 4: Data-to-MC comparison for the multiple-Th final state: (a) the pj distribution of the 
Th candidate in the tt CR; (b) tflj distribution in the Z(— > + jets CR; (c) £x distribution in 
the W + jets CR; and (d) tflj distribution with the requirement | A<p{]2, Vlj) \ < 0.1. The bottom 
panes show the ratio between data and background while the hatched area depicts the total 
uncertainty on the MC. 

6.2.2 Estimate of the Z(— ^ vv) + jets event background to the multiple-Th final state 

The contribution of Z(— > vv) + jets events to the multiple-T^ SR is evaluated by selecting a 
sample of Z(— > + jets events and treating the muons as neutrinos. The sample is collected 
using a a trigger designed to select a muon and a Th- Jet selection criteria similar to those used 
for the SR sample are imposed. In addition, we require two muons passing the criteria outlined 
in Section |3j The control sample has a purity of about 99% as estimated from simulation. The 
tylj distribution for events in this CR is shown in Fig. Hb). The Z(— > vv) + jets background is 
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estimated by interpreting the pj of the pair of muons as #?t- In order to predict the Z(— > vv) 
+ jets rate in the SR, the Z(— > + jets sample is corrected for the ratio of the branching 
fractions R = B(Z — > vv)/B(Z — > jiji), for trigger efficiencies, for the geometric acceptance 
as measured from simulation, and for the reconstruction efficiency e^ eco as measured from data. 
Therefore, CCjj in Eq. Q is given by: 

1 B(Z ^ vv) e^f 8er 



A 2 £ reco2 B(Z ^ w ) £ Trigger 



Since there is no prompt production of a genuine Th in the Z(— > + jets sample, cc x ; = and 
a TT = 0. The Z(— > vv) + jets contribution to the SR is calculated as: 

twCR Trigger 

N |R = ^_W±^ R ^L^Tp (2)/ (8) 

z,^i/i/+jets y^2 £ reco2 Trigger v v ' 

where £^ 88M is the trigger efficiency and £^x 88er the //Th trigger efficiency. The efficiency for 
the #?t > 250 GeV signal selection (e^ 1 ) is determined by calculating the fraction of the observed 
events in the CR that have Vlj > 250 GeV. The muon identification efficiency is measured 
using a "tag-and-probe" method. The probability for a jet to be misidentified as a Th lepton has 
a measured value of / = 0.016 ± 0.002. The estimated Z(— > vv) + jets contribution to the SR is 
determined to be NfL w) = 0-03 ± 0.02. 

6.2.3 Estimate of the Z(— ^ tt) + jets event background to the multiple-T h final state 

The contribution from Z — > tt events is determined with the Z(— > }i}i) + jets CR sample used 
to estimate the background from Z(— > vv) + jets, with the muons treated as Th leptons. The oc xy 
factors are more difficult to estimate for Z — > tt events since there are several ways in which 
Z — > tt events can contribute to the SR: (i) both t^ leptons pass the kinematic acceptance and 
identification criteria; (ii) both Th leptons pass the kinematic acceptance criteria, but only one 
passes the identification criteria; (iii) one t^ fails the kinematic acceptance criteria, while the 
other Th passes both the kinematic acceptance and identification criteria; or (iv) both Th leptons 
fail the kinematic acceptance criteria. The Z(— > tt) + jets contribution to the SR is calculated 
as: 



A 2 T e 2 T 2A 2 T e T (l-e T ) 

/\2„rero2 ' V 1 / 



/\2 P reco2 1 /\2 P reco2 



/\2 P reco2 1 \ X J ^ /\2 P reco2 ' v*v 
■ /L fi c -fi lx \i^]i 



(9) 



where R is given by: 



B(Z^TT)B 2 (T^T h ) gff 

B(Z -4 wi) £ ™s 



e^ T , (10) 



A T is the Th acceptance efficiency, £ T is the Th identification efficiency in this control sample, 
and / = 0.016 ± 0.002. The estimated Z (— > tt) + jets contribution to the SR is determined to 
be N§* , = 0.21 ± 0.13. 

Z(— S'Tt) 
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6.2.4 Estimate of the W + jets event background to the multiple-r h final state 

To select the W + jets CR, the Th isolation requirement, which discriminates between a Th lep- 
ton and other jets, is removed from the SR selection requirements. However, the lack of the Th 
isolation requirement increases the contribution from other backgrounds as most of the back- 
grounds arise because jets are misidentified as a Th lepton. To minimize the contribution from 
tt production, events are required to have no jets identified as a b jet. This requirement reduces 
the contamination from tt events to around 5%. The purity of the W + jets CR is approximately 
65%. Figure Qc) shows the £x distribution, defined as the magnitude of the negative of the 
vector sum of the transverse momentum of all PF objects in the event, for events in the W + jets 
CR. The contributions of QCD multijet, tt, and Z(— > vv) +jets events are subtracted in order 
to determine the number of W + jets events in the CR. The predicted rates for QCD multijet, 
tt, and Z(— > vv) + jets events are determined by extrapolating from their corresponding CRs. 
Since there is no genuine multiple-Th production in W + jets, a TT = 0. According to simulation, 
the fraction of events in the CR with one genuine Th is A T j = 0.149 ± 0.016, while the fraction 
of events without a genuine Th is Ajj = 0.851 ± 0.038. Therefore, u T j in Eq. ^ is given by 
A T j£^° / P(0 b jets), where e 1 ^ is the probability for a Th to pass the isolation requirement and 
P(0 b jets) is the probability to not have any light-quark or gluon jet misidentified as a b jet. 
Similarly, ot;; is given by Ajj/P(0 b jets). The contribution of W + jets events to the SR is then 
calculated as: 



After subtraction 

= piquets) K^ 1 ) + V(2)] • (11) 

The average rate at which jets are misidentified as a Th lepton is measured to be 0.019 ± 0.001. 
The rate /b at which light-quark jets or gluon jets are misidentified as a b jet is used to determine 
P(0 b jets). The estimated W + jets contribution to the SR is determined to be N^_- ets = 5.20 ± 
0.63. 

6.2.5 Estimate of the QCD multijet event background to the multiple-r h final state 

QCD multijet events contribute to the multiple-Th SR when mismeasurements of jet energies 
lead to large values of Vlj and when jets are misidentified as Th candidates. By removing the 
Th isolation criteria and inverting the |A</>(j2, #?t) | requirement, a QCD CR sample with about 
99% purity is obtained. Figure |4]jd) shows the expected and observed Iflj distributions for this 
sample. A scale factor is obtained from this CR and used to correct the signal prediction for 
QCD multijet events in simulation. The estimated contribution to the SR from QCD multijet 
events is determined to be N^ D = 0.02 ± 0.02. 



7 Systematic uncertainties 

Systematic uncertainties are taken into account for both signal and background events and are 
described separately. Both the signal and background are affected by the systematic uncertainty 
in the identification of the Th candidate. The systematic uncertainty for Th identification is 
obtained using a Z — > tt enhanced region and by correcting this cross section by that measured 
for Z — > ee and Z — > fi}i events. This uncertainty is validated on a control sample of Z — > tt 
events. The level of agreement between data and simulation is found to be at the level of 
7%. Further validation of the performance of Th identification in a SUSY-like environment is 
performed by selecting a W(— > tv — > t^vv) + jets CR with large hadronic activity (Hj) and large 
transverse momentum imbalance {Vij). The level of agreement between the predicted rate for 
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W — > zv — > TfrW events and the observed number of events is within 7% and is determined as 
a function of Hj and t/tj. 

7.1 Systematic uncertainties on background events 

The principal sources of systematic uncertainty on the background predictions arise from the 
correction factors, the finite number of events in the CRs, the measured rates at which jets are 
misidentified as a Th lepton, and the level of agreement between the observed and predicted 
numbers of events in CRs. 

The contributions to the uncertainties on the correction factors are different for each back- 
ground category. The dominant effect is due to the uncertainty in the Th identification efficiency. 
In the multiple-Th final state, uncertainties in the jet-energy scale (JES) |32| and the Th-energy 
scale (TES) II33II are used to evaluate how changes in Hj, Iflj, and jet kinematics affect the cor- 
rection factors. The systematic uncertainty on the correction factors due to the JES and TES is 
at most ~3% for all backgrounds. Smaller contributions to the uncertainties in the correction 
factors arise from the muon reconstruction and isolation efficiency (<1%), the uncertainty in 
the branching fractions (<Cl%), and the uncertainties in trigger efficiency (1%). 

The systematic uncertainties on the measured rates for jet misidentification as a Th lepton are 
dominated by the size of the jet sample used to measure these rates and range from 2% for the 
single-Th final state to 5.6-10% for the multiple-Th final state. The level of agreement between 
the observed and predicted number of events in MC studies of the CRs is used to assign an 
additional systematic uncertainty and ranges from 2% for the single-Th final state to 3% for the 
multiple-Th final state. Finally, the systematic uncertainty arising from statistical uncertainties 
on the number of events in the CRs ranges from 2-5% for the multiple-Th final state to 3-10% 
for the single- Th final state. 

7.2 Systematic uncertainties on signal events 

The main sources of systematic uncertainties in the SR are due to trigger efficiencies, identifi- 
cation efficiencies, the energy, and momentum scales, the luminosity measurement and PDFs. 
The uncertainty on the luminosity measurement is 2.2% [34|. Systematic uncertainties on the 
#?t triggers (2.5%) are measured using a sample in which around 99% of the events are tt events, 
which have a similar topology to events in the SR samples. The systematic uncertainties on the 
TES and JES (3.0%) yield an uncertainty on the signal acceptance of 2.3%. The uncertainty on 
the yij scale depends on the uncertainty of the JES (2-5% depending on the t] and pj values of 
the jet) and on the unclustered energy scale (10%). Unclustered energy is defined as the energy 
found "outside" any reconstructed lepton or jet with pj > 10 GeV. The unclustered energy 
scale uncertainty has a negligible systematic uncertainty on the signal acceptance. The system- 
atic uncertainty due to imprecise knowledge of the PDFs (11%) is determined by comparing 
the CTEQ6.6LE3, MSTW 2008 NLO HO, and NNPDF2.1 |37| PDFs with the default PDF HQ. 
The systematic uncertainty due to the imprecise modeling of the initial-state and final-state 
radiation [39| is negligible (<Cl%). The systematic uncertainties associated with event pileup 
are also negligible. Uncertainties on the theoretical cross sections are evaluated by varying the 
PDFs and by changing the scales by a factor of two |25l - l3"0fl . 

8 Results 

For the single-Th final state, the number of background events containing a genuine Th, as well 
as the number of background events containing a misidentified Th, are estimated with data. 
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The results for the baseline and the full selection are listed in Table [3] Figure [5] shows the Hj° 
and yij distributions of data and the different background predictions. The observed number 
of events in data is in agreement with the SM predictions. 

Table 3: Number of data and estimated background events with statistical and systematic un- 
certainties, respectively, in the single- Th final state. 



Process 


Baseline 


Signal Region 


Fake— Th 


67±2±19 


3.4 ± 0.4 ± 1.0 


Real— Th 


367 ± 10 ± 27 


25.9 ± 2.5 ± 2.3 


Estimated £ SM 


434 ± 10 ± 33 


29.3 ±2.6 ±2.5 


Data 


444 


28 




Figure 5: Distributions of (a) H| , and (b) tflj for the single- Th final state. The points with errors 
represent data that satisfy the baseline selection while the filled green (light) and filled blue 
(dark) areas shows the predicted backgrounds due to events containing a genuine Th and a 
misidentified Th, respectively. The hatched area shows the total uncertainty on the prediction. 

The largest sources of background for the multiple-Th final state are from tt and W + jets events. 
A counting experiment is performed and the background predictions from data are compared 
with the observed number of events. Table [4] lists these background predictions and the ob- 
served number of events in the SR. Figure [6] shows the Hj° as well as the M e ff distributions in 
the SR. The background distributions in Fig. |6]are taken from simulation and normalized over 
the full spectrum. The estimated number of events due to the SM background processes is in 
agreement with the number of observed events in the SR. 



9 Limits on new physics 

The observed numbers of events in the single-Th and multiple-Th final states do not reveal any 
evidence of physics beyond the standard model. Exclusion limits are set using the CL S EUI 
criterion in the context of the the CMSSM |4T|. The CMSSM parameter space with tan j6 = 40, 
Aq = — 500 GeV, ]i > 0, and M t = 173.2 GeV is chosen as a possible scenario with a light f 
and a value of AM < 20 GeV. The excluded regions are shown for the single-Th and multiple- 
Th final states in Figs. [7ja) and |7jb), respectively. The limits are set using a simple counting 
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Table 4: Number of data and estimated background events with statistical and systematic un- 
certainties, respectively, in the multiple-Th final state. 

Process Signal Region 



QCD multijet events 0.02 ± 0.02 ± 0.17 

W+jets 5.20 ± 0.63 ± 0.62 

tt 2.03 ± 0.36 ± 0.34 

Z(^TT) + jets 0.21 ±0.13 ±0.17 

Z(-M/v)+jets 0.03 ± 0.02 ± 0.50 

Estimated £ SM 7.49 ± 0.74 ± 0.90 

Data 9 
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Figure 6: Stacked distributions of (a) H|°, and (b) M e g in the SR for the multiple-Th final state. 
The background distributions are taken from MC events that are normalized to the predictions 
based on data over the full region. The shapes obtained from MC simulation are used for 
illustrative purposes only. 



experiment. Systematic uncertainties are treated as nuisance parameters and marginalized, 
and the signal contamination is taken into account. In the CLs method, both the background- 
only as well as the signal + background hypothesis are used to derive the confidence levels CL S 
and the resulting limits and the uncertainty bands on the exclusion contours. In the case of very 
small values of AM(~5GeV), the lower-energy cannot be effectively detected and only the 
energetic from the decay of the neutralino can be observed. The search for new physics with 
a single t lepton has a better sensitivity in this case. The single- and multiple-Th topologies 
thus have complementary sensitivity and together provide coverage for models with a wide 
range of AM values. 

Using the limits set by the single- analysis, a common gaugino mass HI1/2 of <495GeV is 
excluded at 95% Confidence Level (CL) for a common scalar mass niQ of <400 GeV. For the 
multiple-Th analysis, m\/2 < 495 GeV is excluded at 95% CL for mo = 400 GeV. A gluino with 
mass <1.15TeV is excluded at 95% CL for mo < 500 GeV. It can be noted that the single-Th 
analysis shows better sensitivity for small values of AM, which is near the boundary of T = 
LSP 

The results for the multiple-Th final states are also interpreted in the context of SMS |13|. The tt 
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Figure 7: (a) Exclusion limits in the CMSSM plane at tan = 40 for the single-ih final state. The 
observed and expected limits are shown by the solid red and blue lines, respectively. The blue 
band represents the expected uncertainties while the red dotted lines represent the theoretical 
uncertainties for events that satisfy the full selection, (b) Exclusion limit in the CMSSM plane at 
tan /3 = 40 for the multiple-Th final state. The solid red line denotes the experimental limit while 
the dotted red lines represent the uncertainty on the experimental limit due to uncertainties on 
the theoretical cross sections. The blue band represents the expected uncertainties. 
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SMS scenario (T3tauh) is studied where gluinos are produced in pairs and subsequently decay 
to t lepton pairs and an LSP via a neutralino (g — > qq^; — > tt — > tt^). The diagram for 
the T3tauh model is given in Fig. [HJ A gluino mass of <740 GeV is excluded at 95% CL for LSP 
masses up to 290 GeV (here, the mass of is the average of the masses of the gluino and the 
LSP). Figure [9jja) shows the 95% CL exclusion region obtained for T3tauh. The limits on the 
mass of the gluino and LSP are shown with a solid black line. 



q 




Figure 8: Diagram for the T3tauh SMS model. 
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Figure 9: Exclusion limits for the multiple- final state: (a) 95% CL cross section upper limits 
for the T3tauh model, where the solid red line represents the limits on the mass of the gluino 
and the LSP; (b) 95% CL cross section upper limits as a function of gluino mass in the GMSB 
scenario. 



In the simplified Gauge Mediated Symmetry Breaking (GMSB) scenario, the t is the NLSP 
and decays to a T lepton and a gravitino G, with a mass of the order of ~ keV lj42Tl44| — > 
tt —¥ ttG). The topology for this simplified GMSB scenario is similar to that of T3tauh except 
for the assumption that both the gluinos decay to T-lepton pairs with a branching fraction of 
100%. Therefore, the results are also interpreted in the simplified GMSB scenario using the 
T3tauh scenario. The signal acceptance is corrected to account for the final state containing up 
to four r leptons. A gluino with mass <860 GeV is excluded at 95% CL. Figure |9pb) shows the 
exclusion limits for the simplified GMSB scenario as a function of the gluino mass. The ATLAS 
collaboration has also reported a result on a search for one or more hadronically decaying tau 
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leptons, highly energetic jets, and a large transverse momentum imbalance probing minimal 
GMSB models g5]|. 

Since the SMS topologies considered in this paper are characterized by two x leptons in the 
final state, we do not present SMS limits for the single-T^ final state. 

10 Summary 

A search for physics beyond the standard model with one or more hadronically decaying r 
leptons, highly energetic jets, and large transverse momentum imbalance in the final state is 
presented. The data sample corresponds to an integrated luminosity of 4.98 ± 0.11 fb _1 of pp 
collisions at y/s = 7 TeV collected with the CMS detector. The final number of events selected 
in data is consistent with the predictions for standard model processes. We set upper limits on 
the cross sections for the CMSSM, GMSB, and SMS scenarios. Within the CMSSM framework 
at tan /3 = 40, a gaugino mass rn\/2 < 495 GeV is excluded at 95% CL for scalar masses mo < 
400 GeV. This result sets a lower limit on the mass of the gluino at 1.15 TeV with 95% CL in this 
region. In the multiple-Th final state, a gluino with a mass less than 740 GeV is excluded for 
the T3tauh simplified model while a gluino with a mass less than 860 GeV is excluded for the 
simplified GMSB scenario at 95% CL. 
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